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Abstract. We review the application and successes of a phase-space coalescence 
plus fragmentation model that has been applied for hadronization at RHIC. The 
physical concept is discussed together with the practical implementation. The ro- 
bustness of main predictions is reviewed together with several open issues like 
relevance of three dimensional calculation, finite width of the wave functions, ef- 
fects of quark masses, energy-entropy conservation, space-momentum correlation. 
Eventually the relevance of coalescence also for the study of the microscopic in- 
teraction of heavy quarks is highlighted. 

1 Introduction 

The goal of ultra-relativistic heavy-ion collisions is to create a new state of nuclear matter with 
quarks and gluons as degrees of freedom. Such a state of matter is predicted to be one of the 
Quantum Chromodynamics (QCD) phases, usually referred to as quark-gluon plasma (QGP), 
and it is expected to occur at energy densities above tc — 0.7 GeV/fm"^ and temperatures above 
Tc ^ 180 MeV [1]. These conditions are believed to be reached at the Relativistic Heavy Ion 
Colhder (RHIC) as confirmed by the success of ideal hydrodynamics [2l3j (at lowpr <l-5 GeV) 
and the observed strong nuclear suppression Raa in agreement with jet-quenching predictions 
[3]. In fact both entail an energy density e > 10 GeV/fm'^ ^ e^- 

Such a hot and dense medium is therefore expected to be in a deconfined state, i.e. a state 
where quarks and gluons act directly as degrees of freedom. The success of hydrodynamics 
based on a QGP equation state in reproducing the observed large elliptic flow ^2 3 , together 
with the failure of hadronic transport model [5] , constitutes a first indirect proof of the partonic 
phase. However in the search for signatures of the QGP and in the study of its properties an 
important role should be played by the understanding of hadronization dynamics, as persistently 
stressed by J. Zimanyi [5]. However because of its intrinsic non-perturbative nature there has 
been a general skeptical view on the real possibility to have hints of quark degrees of freedom. 
Surprisingly the wide availibility of exclusive data from RHIC experiments and the power of 
the phenomenological approach have made possible to spot signatures of hadronization via 
coalescence of massive quarks I7|8|9|10|ll|12j . 

Coalescence is a quite general concept that has been applied to nuclear physics for light clus- 
ter production for more than two decades |14j . First application to ultra-relativistic heavy- ion 
collisions is due again to J. Zimanyi and his group in Budapest who implemented an ALge- 
braic COalescence Rehadronization (ALCOR) model that studies the chemical composition of 
the hadronic matter produced in heavy ion collisions at SPS and RHIC T5]. However it was 
only in 2002-03 that some puzzling behavior of the experimental data in Au+Au collisions 
at RHIC showed features specific of valence quark coalescence: enhancement of baryon/meson 
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ratio at intermediate pr, absence of Raa suppression for baryons, and the scaling of elliptic 
flow according to the costituent quark number |8|9|ll|12|13j . 

In this talk, we will review the general idea of coalescence and its phase-space implemen- 
tation for the hadronization. Then we will show the comparison with the experimental data 
for the particle spectra, their ratios and elliptic flows. Finally we will discuss the robustness 
of the basic feautures of a naive coalescence approach together with some open issues. Quark 
coalescence was initially applied to light quarks, but I will also draw the attention to its role 
in the investigation of the in-medium interaction of heavy quarks and quarkonia in the QGP 
|10I16I17| . 

2 Coalescence and Fragmentation 

A lot of efforts has been done to work on the hadronization process, even if a fully satisfying 
theoretical description is still missing. Because of its highly non perturbative nature, hadroniza- 
tion is treated through different phcnomenological approaches that have shown their validity 
in different momentum scale regions: string fragmentation, indipendent fragmentation, dual 
parton model, coalescence. We will discuss mainly a coalescence plus fragmentation model that 
has the advantage to describe in an economical way different features of the hadronization and 
the collective behavior of the hadronic matter created at RHIC. 

We briefly remind the main features of indipendent fragmentation (IF) used as standard 
hadronization mechanism for enough large pT- IF relies on the QCD factorization theorem that 
allows to write the hadron production as a convolution /h(pt) = fpiPT / z) (8) Dp^niz, fx), 
where fp is the parton distribution in the system, Dp^h is the fragmentation function, with 
z fraction of the momentum of the parton carried by the hadron and fi ^ the pertuba- 
tive scale. In e+e^ and hh' collisions, fp can be determined by the pQCD cross section folded 
with the parton distribution functions (PDF's), while in heavy-ion collisions they are calcu- 
lated as sovrapposition of pp collision plus nuclear shadowing and, moreover, radiative energy 
loss due to in- medium gluon-radiation In such an approach non-pertubative features are 
encoded in the fragmentation function Dp_>/f(z, /i), that gives the probability to produce from 
a parton p a hadron with momentum ph — zpx- Such a function is considered to be univer- 
sal and is extracted from the available data in e~^e~ collisions and then applied to study also 
hadron-hadron collisions. This approach has been successful in reproducing the 7r° transverse 
momentum distribution down to pt ~ 2 GeV [22] for pp collisions at RHIC energy. 

Despite the success of IF for high energy collisions and hadrons at large pr, at lower energy 
or large rapidity there have been evidences of hadronization mechanism that proceeds through a 
coalescence mechanism. For example the D~ / ratio in hadron- nucleus reaction at Fermilab 
[20|21j or the particle production in the fragmentation region (large rapidity) that cannot 
be explained by IF while coalescence seems a more suitable approach }19j . In such physical 
conditions it seems necessary to postulate that hadrons come from a convolution of two parton 
distribution folded by a wave function fnip) = /p(pi) '9 fpiP2) ®'^m{pi— P2), where pi and p2 
are the momenta of coalescing partons. One naively expects that when the phase space is enough 
dense and the process involves a low virtuality, the production of quarks from vacuum is not 
likely and the recombination of quarks becomes the dominant mechanism of hadronization. In 
this perspective HIC's provide the ideal environment much denser than hp collisions in a wider 
kinematical region (respect to hh' collisions), for which is natural to think that recombination 
processes may play a dominant role, at least at not much high px- 

On the other hand it is quite natural to think that nature has a smooth transition between 
the two processes as a function of phase space and virtuality, hence IF or, at lower pT^ string 
fragmentation could always be present together with coalescence. 

3 Phase -Space Formulation and Implementation 

We describe our implementation of coalescence here, while the fragmentation part is described 
in the next subsection. Our approach is based on the Wigner formalism [M] that allows a more 
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direct connection with the dynamical phase-space description of heavy- ion collisions (HIC's). In 
this formalism the transverse momentum spectrum of hadrons that consist of n (anti-) quarks is 
given by the overlap between the hadron wave function and the n quark phase-space distribution 
function fq{xi,pi): 



-^^=9H J n (2l)3£:- ^' ■daJq{x,,p,)fH{xi..Xn;pi--Pn)S'-'^^ (pt-^PtA (1) 

where da denotes an element of a space-like hypersurface, fH{xi..Xn',Pi--Pn) is the Wigner 
distribution function of the hadron, gn is the probability of forming from n colored quarks a 
color neutral object with the spin of the hadron considered. In Eq.lU]) it is already assumed 
that the n quark phase space distribution is approximated by the product of the single quark 
distribution function 

n 

fqixi..Xn;Pl-Pn) ^'[[f{x^,P^) (2) 

1=1 

and therefore no quark-quark correlations are included. The formalism can be formally extended 
to include such correlations 27\, but then it is necessary a dynamical study of the in medium 
correlations that asks for transport approaches. 

As light hadrons wave function we have used a sphere in both space and momentum, with 
radii Ar and Ap, respectively, which in the Wigner formalism are related by Ar ■ Ap = 1. A good 
description of pion, kaon, proton, antiproton spectra can be obtained with a radius parameter 
Ap = 0.24 GeV for mesons and 0.36 GeV for baryons, which in terms of mean square radius 
corresponds to take a slightly larger radius for baryons respect to mesons. The multidimesional 
integral Eq.(IT]) is evaluated in the full 6D phase space by the Monte Carlo method via test 
particle method ^9j. 



3.1 Bulk properties 

For partons in the quark-gluon plasma we take a thermal distribution for transverse momenta 
up to po = 2 GeV in agreement with the hydrodynamical behavior observed at low px- For 
their longitudinal momentum distribution we assume boost-invariance, i.e. a uniform rapidity 
distribution in the rapidity range y G (—0.5,-1-0.5). To take into account collective flow of 
quark-gluon plasma, these partons are boosted by a flow velocity vt = Po{rT/R), depending 
on their transverse radial position rx, R being the transverse size of the quark-gluon plasma 
at hadronization, and Po is the maximum collective flow velocity of the quark-gluon plasma. 
Therefore for light quarks and antiquarks transverse momentum spectra are given by 

dNq^q _ gq.qTniT / 7T(mT - PT ' vt T Mg) 

d^VTd^PT ~ (2^)3 ^ T 

where gq = gq = 6 is the spin-color degeneracy of light quarks and antiquarks, and the mi- 
nus/plus signs of chemical potentials are for quarks and antiquarks, respectively. The slope 
parameter T is taken to be T = 170 MeV, consistent with the phase transition temperature 
from lattice QCD calculations p. 

The masses of thermal quarks are taken to be those of constituent quarks, m„ ^ = 300 MeV, 
fUs = 475 MeV, rUc = 1.4 GeV and mi, = 4.8 GeV (effect of masses on the typical observed 
scaling is discussed in Section[4]). We notice that in such an approach non perturbative effect are 
also encoded in the quark masses, in fact a large part of the interaction can be accounted for by 
quark thermal masses |23I24I25] : however the relation of this masses to the masses associated 
to scalar condensates remains still an open challenging theorethical question. For the quark 
chemical potential fj,q, we use a value of /ig = 10 MeV to give a light antiquark to quark ratio of 
0.89, which would then lead to an antiproton to proton ratio of about (0.89)^ = 0.7, consistent 
with the observed ratio at midrapidity in heavy ion collisions at RHIC. 
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Above po, partons are from the quenched pQCD minijets [18], and their masses are those 
of current quarks (which schematically resembles the p— dependence of the constituent quark 
masses). Both soft thermal and hard minijet partons are assumed to be distributed uniformly in 
a fireball of volume of 950 fm'^ (which implies a transverse radius R = 8.1 fm and longitudinal 
length of r = 4.4 fm for a Bjorken plus transverse expansion). Such a volume is fixed to 
reproduce the measured transverse energy of 750 GeV |34| , together with radial flow parameter 
(3o = 0.5 which is consistent with both experimental data and hydrodynamical calculations [3j. 
In other words the reconstruction of the fireball tells us that the system used for coalescence 
hadronizes at an energy density of 0.8 GeV/fm"^, which is essentially the energy density at 
which the phase transition is expected to occur from lattice QCD calculations T. In addition 
the entropy of such a system is dS/dy = 4800 in agreement with the value inferred from 
experimental data by S. Pratt and S. Pal |45| . Therefore even if employed mainly at intermediate 
Pt coalescence model relies on a bulk distribution that is consistent with what can be inferred 
from hydrodynamics and experimental data. 

3.2 Minijet Distribution 

The transverse momentum distribution for px > 2 GeV is taken to be that of minijet partons 
in the midrapidity. In particular we use the one that can be obtained from an improved pertur- 
bative QCD calculation [3Dj. It is given by dNj^t/d'^pT = ^/'^tot^'^^'^ict/d'^pT in terms of cTtot^" 
corresponding to the total cross section at central 10% of the collisions and the jet production 
cross section from nucleus-nucleus collisions, 



In the above, iA(r) is the thickness function of Au at transverse radius r. The parton distribution 
function in a nucleon in the nucleus Au is denoted by fa/ Au{x, Q^) plus a transverse smearing 
(7(kT). The cross section da°'^ /dt is the parton scattering cross section. Kinematic details and 
a systematic analysis of pp collisions can be found in Ref . [30] . Using the GRV94 LO result 
for the PDF [3T] and the KKP fragmentation function from Ref. [35], measured data in the 
reaction pp tt^X at ^/s = 200 GeV can be reproduced with Q = 0.75pT and {k\) = 2 GeV^. 

In heavy ion collisions at RHIC, minijet partons are expected to lose energy by radiating soft 
partons as they traverse through the quark-gluon plasma [Ij. This effect is taken into account 
by lowering their transverse momenta by the energy loss AE, which depends on both the parton 
energy E and an effective opacity parameter L/\ according to the GLV model [T^. An effective 
opacity L/X — 3.5 is used, as extracted from a fit to the spectrum of high transverse momentum 
pions measured at RHIC [18] 0. 

We note that our approach is indeed quite schematic in separating the thermal spectrum 
from the one of minijets. In reality there is only one spectrum with possible correlations also at 
Pt < Po- Our simplificaton does not affect the one-body observables described in the following, 
but can be essential for studying two or three particle angular correlations [76|75j . 

4 Coalescence for light quarks 

In the coalescence model, hadrons are formed from quarks that are close in phase space if a 
hadron wave function with a small width is considered. As a result, baryons with momentum 

^ Such an assumption is in principle too drastic for minijets, however it does not affect our results 
because jet-jet coalescence is never taken into account. 

^ In this calculation is not included the QCD analog of the Ter-Mikayelian effect. When included a 
L/A = 5 has to be used to have the same amount of quenching [61] . 




X fa/Au{Xa,Q'^)fb/Au{xb, Q'^)-6{s + t + u) 



da' 



(4) 



dt 
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Fig. 1. Left: Pion and antiproton transverse momentum spectra from Au+Au collisions at y^s — 200 
AGeV, calculation from coalescence plus fragmentation of minijets is shown by solid lines; experimental 
data [35] for 7r° are shown by filled circles, and for p by squares. Right: p/n and K^/A ratios, model 
by solid lines and experimental data [33138) by filled squares. 

Pt are produced from quarks with momenta ~ PtI'^i while mesons with same momentum are 
from quarks with momenta ^ pt/'2,. Since the transverse momentum spectra of quarks decrease 
with Pt, production of high momentum baryons from quark coalescence is favored respect to 
the fragmentation where baryons are penaUzed with respect to mesons as more quarks are 
needed from the vacuum. In the independent fragmentation the estimated value of p/tt ratio is 
of ~ 0.2 0.3 in approximate agreement with the ratio measured in pp collisions. For Au+Au 
collisions results for p and tt" spectra are shown (solid lines) together with the data from 
PHENIX in Figlll^left). The resulting p/ir and the K^^/A ratios based on Eq.dT]) are shown 
in Figd] (right) together with data points [26 38 . Calculations shown for tt, K include the 
constribution from indipendent fragmentation that is essential to reproduce the meson spectra 
already dXpT ^ "i GeV (especially for pious). Although different models have been used, they all 
lead to enhanced p/tt ratio and K^/A ratio of similar magnitude |7|llj . As shown in Fig[T](left), 
our approach [8*9] , which includes resonance decays and avoids the coUinear approximation by 
using the Monte Carlo method to evaluate the multi-dimensional coalescence integral, gives a 
good description of spectra, p/n and K'^/ A ratios also at < 2 GeV, even if in this momentum 
region the assumptions used are less under control (for a more extended discussion see next 
Section). At lower energy, Au+Au at 62 AGeV, a direct extrapolation with the same model 
predicts [66] an enhancement of thep/7r ratio and a decrease of p/tt respect to 200 AGeV. Such 
a trend is again in reasonable agreement with experimental data [28129] . 

In HICs the anysotropy of particle momentum distributions with respect to the azimuthal 
angle offers the possibility to get information on the dynamics of the collisions and the properties 
of the created hot and dense matter i40j . At RHIC energies it has been shown to be an important 
probe for the equation of state of the quark-gluon plasma[3] and of the parton cross section 
|36I37| , and for hadrons at high transverse momentum it is a probe of the initial energy density 
that causes the jet quenching Moreover it has been shown that it is built-in very early in 
the dynamical evolution being thus a probe of the interaction during the partonic stage [36] , a 
picture recently confirmed by the measurement of meson V2 |51| . 

Thanks to coalescence models it was realized that elliptic flow provides also a way to under- 
stand the hadronization mechanism in itself. In fact if hadronization goes through coalescence 
of constituent quarks the anisotropy at partonic level propagates at hadronic level according 
to[l2] : 

V2,m{pt) ~ 2w2,g(pT/2) , V2,b{pt) ~ 3u2,9(pt/3). (5) 

This effect together with the shape of rise and saturation at parton level predicted by the 
parton cascade [37] for the partonic flow allows to see the charateristic feature of coalescence 
that is a larger elliptic flow for baryons, by a factor 3/2, and a similar shift of the pt at which 



6 



ditor 




Fig. 2. Left: quark number scaled elliptic flows of tt and p in Au+Au at 200 AGeV. Lines are from 
the coalescence model, while symbols are data from PHENIX [42j. Right; (top) scaled elliptic ow for 
identified hadrons. (bootom) ratio of difference between quark number scaled baryon V2 and quark 
number scaled meson V2 divided by the sum: {B — M)/{B + M). Model predictions are also shown on 
the lower panel. Figure is taken from Ref. [65) . 



the V2 reaches the maximum value. These features are clearly seen in the experimental data, 
as anticipated by Voloshin [13] , see FiglS] We notice that the coalescence mechanism tanslates 
the hydrodynamical behavior into higher px at hadronic level with the effect being larger for 
baryons. This is seen also from the agreement of hydrodynamics V2{pt) for mesons up to 1.5 
GeV and for baryons up tp ^ 2.5 GeV . As a general remark we like to point out that coalescence 
is the only known hadronization mechanism that can lead to an enhancement of the partonic 
elliptic flow. Someone may think that hydrodynamics gives the larger V2, but at fixed elliptic 
flow in the partonic stage coalescence can still enahance the elliptic flow. However, Eq.(l5]) is 
valid under the assumption of a uniform phase space density |46| and the approximation of a one 
dimensional space in which only coUinear and equal momentum partons reconibine. Relaxation 
of these assumptions can in principle break "quark number scaling" of the elliptic flow. For a 
discussion about effects beyond the simplest (naive) coalescence, see the next Section. 

In Fig. [5] (left) we show the scaled elliptic flow for pions and protons obtained with our model 
once the quark V2q is fitted to reproduce K elliptic flow, then elliptic flow of p, p, K, A is well 
reproduced |9|68j . The calculations include the effect of resonance decays that is another possible 
source of scaling breaking |50| , as dicussed more in detail in the next Section. The contribution to 
hadrons elliptic flow from minijet fragmentation is not included in this calculation. Its inclusion 
would lead to a universal hadron elliptic flow at momentum above 6 GeV (llj . 

In Fig. [2] (right), taken from Ref. [65] . the scaling of V2 for various identifled hadrons is shown 
together with the relative baryon to meson {{B — M) / {B + M)) deviation from the naive scaling 
Eq.(l5]). The circles are the experimental data, the dotted line is the prediction of our model 
|50j , and the solid line is the effect coming from higher Fock state in the hadron wave function 
discussed in Ref. [78\ . The agreement of our model with the data comes from a combination of 
the finite width of the hadron wave function and the effects of resonance decays. 

In FigO (left) it is shown the effect of quark mass on the scaling of V2 at fixed wave function 
width (Z\p = 0.24 GeV). In panel (a) it is shown an ideal case in which the quark V2q{pT) is 
flat, it is evident a strong violation of the scaling at low pr- The source of such a violation of 
the scaling is the full 6D phase space that allows coalescence of non-collinear momentum at low 
Pt where the boost effect is small. In such a case the effect is nearly independent on the quark 
mass. In the lower panel (b) a more realistic case is shown, i.e. a V2{pt) that rise and saturates 
at V2 = 0.1. In such a case the violation of the scaling at low pT is drastically reduced, while 
at higher pT it is significantly mass dependent. We clearly see a violation of ~ 25 — 30% if a 
mass rUq = 0.03 GeV is used. This is due to the fact that when the mass is small, even if the 
wave function width is small in the rest frame of the hadron, particles with large difference in 
the relative momenta can still coalesce due to boost effects. In such a case particles at higher 
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Pt (w2 ~ 0.1) can coalesce with particles at small pT (w2 ^ 0) and this effect is of course larger 
for smaller nnark ma,ss. 




p^/n [GeV] (GeV) 



Fig. 3. Left: Scaled elliptic flows, the solid black line is the quark V2q-, the dashed and dot-dashed 
line are the V2 of mesons from coalesccence for two different values of the quark mass ruq, see text for 
details. Right: V4 in Au+Au collisions at 200 AGeV for charged (solid line), pion (dashed line) and 
proton (dot-dashed line) evaluate from coalescence; circles are data taken from Ref . [53) . 



It will be important to see experimentally where the scaling breaks down as a function of the 
beam energy but also as a function of the rapidity. In principle one would expect that at very 
forward (backward) rapidity one reaches conditions that are comparable to the midrapidity 
one at lower energy, therefore it could be that a QGP phase is no longer dominant at large 
rapidity. Indeed we have seen that in the AMPT model the magnitude of vi and W2 at y > 3 
are consistent with a pure hadronic matter dynamics [43147] . 

The study of azimuthal anisotropy has been mainly limited to the elliptic flow (second 
harmonic of a Fourier expansion in the azimuthal angle) but a sizeable amount of the fourth 
order momentum anisotropy has been predicted '54] and confirmed experimentally [52 . This 
provides the possibility to further test the coalescence model, in fact also for the va precise 
relations between mesonic and baryon V4 are expected 39j, as for example: 

-^(3pt)« :r-^(2pT) + ^ (6) 

Finite momentum spread and resonance decay do not change significantly such a prediction, 
especially at pt > 1 GeV [68139] . For the details on the derivation of this and other scaling 
relations also for odd harmonics we remand to Ref. [55.. In Figl3| (left) predictions for pion and 
proton Vi are shown, the calculation is done fitting to the charged hadron W4. Recent data 
[56] again confirm the trend predicted by coalescence [39 68 . 



5 Beyond the simplest implementation 

The RHIC program has provided a remarkable evidence that coalescence of massive quarks 
supplies a simple and succesful model for hadronization from a deconfined plasma. Nonetheless 
there are several aspects that are still problematic and act as a stimolous to a deeper and 
advanced formulation of the coalescence process in the context of hadronization. Some of these 
problems are more serious if hadronization via coalescence is applied to bulk production of 
hadrons. Nevertheless, once coalescence is recognized to be the dominant mechanism at inter- 
mediate Pt, it is appealing to investigate its behavior at low momenta where the phase space is 
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Fig. 4. Left: Scaled V2 for (a) pions coming from different resonane decays and (b) for K and p from 
K* and A respectively. Dotted line show the underlying quark V2q- Right: scaled V2 for quarks (solid 
line) and for mesons for different wave function width Ap. 



denser and coalescence probability increases. Till now we have simply extended the approach 
at quite low momenta noticing that no striking contradiction with the experimental data are 
observed which validates also the good agreement at intermediate pr- Even if one should be 
aware of the accompanying problems, that we discuss in the following, also the success of AL- 
COR in reproducing the total yields [T3] has to be considered. Furthermore the reasonable 
description of the experimental data on particle spectra also at low px in our approach is found 
also in a more recent based on ALCOR/MICOR approach [44]. We want also to point out 
that an important probe for the formation of QGP and the hadronization mechanism is given 
by the charge fluctuations '72]. First data fT? agrees better with predictions from coalescence 
[73j . Furthermore charge fluctuation of the bulk are reproduced by a coalescence hadronization 
mechanism if the number of quarks and antiquarks is of the order of dN/dy = 1300 [5T] (for 
most central collisions) which is similar to the number used in the model discussed here at 
intermediate pr and to what is extracted by ALCOR model that is instead dedicated to the 
study of the multiplicities [80]. This shows again an internal consistency of coalescence models. 

The simplest approach to hadronization by coalescence considers a uniform distribution of 
particles that combine if have the same pr forming directly the stable hadrons. Of course this 
seems too naive for a number of reasons that I try to list here together with some investigations 
beyond this simple picture that have been or could be done in the next furture: 

1. Resonances - The baryon/meson anomaly and moreover the scaling of elliptic flow are 
directly applicable to stable hadrons, on the other hand it is well-known that especially for 
pions there is a large feed-down from resonance decays. In the model we have presented here 
and published originally in Ref. [9| resonances were already included. As discussed more in 
detail in Ref. [50] , at intermediate pt their role is quite reduced and therefore their discard 
by other groups is justified when looking at — 3 GeV. However our calculation shows 
that when included both pT spectra and V2{pt) have a better agreement with data at lower 
Pt- In particular, we have studied how the elliptic flows of pions and other stable hadrons 
are affected by decays of resonances, such as p ^ 2 tt, a; ^ 3 tt, K* ^ Kit, and Z\ — > p7r 
[50] , see FigHJ It turns out that particles like p. A, and K from resonance decays have 
elliptic flows that are very similar to the directly produced ones. Therefore, the inclusion of 
resonances does not destroy the coalescence scaling of these stable hadrons. On the other 
hand, pions from the decay of K* , A show a significant enhancement of their elliptic flow 
at Pt < 2 GeV. The breaking of coalescence scaling due to resonance decays together with 
that due to finite quark momentum spread (see next point) lead to a better agreement with 
available data as shown in Fig|2] where the decay of resonances is included 
However pions are not expected to be well described, but one has also to consider that at 
intemediate pt the mass mismatch is less relevant and at low momenta most of the pions 
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Fig. 5. Left: Evolution of the total entropy in \y\ < 0.5 associated to u,d matter as a function of time, 
for different options about the expansion and the decay of hadrons formed by coalescence. 

are coming from resonance decays. This can be the reason behind the approximate good 
description of pions in our model. [9150] . 

2. Wave function - Our formahsm described in Section [3] contains the full 3D phase space 
and a wave function with a finite width Ap of the order of the Fermi momentum [50j . 
The effect Ap is shown in Fig. [4] (right) where we can see a significant increase of the 
violation of the scaling between partons and mesons as a function of the width of the wave 
function. However with the same width used to reproduce the px spectra the inclusion of 
a full 3D phase space and wave function does not destroy but even improve the agreement 
with experimental data, see FigOand the relative discussion 0. An extreme ansatz for the 
wave function could be considered, but if the wave function is flat in term of the relative 
momentum, the naive formulas Eq.([5]) (right) are strongly modified leading to a reduced 
difference between baryons and mesons |48j not observed experimentally. 

3. Energy and Entropy Conservation - Another issue related to Eq.([T]) is the energy conser- 
vation. In the present approach on shell massive quarks are recombined into one particle 
respecting momentum conservation, but not energy conservation. Of course this is a draw- 
back of the simplest implementation while recombination processes always happen in nature 
involving a third particle or off-shell recombining particles. 

In our model pions come mainly from p decay hence the energy violation is less serious 
even at low px, at the same time both the px distribution and the elliptic flow get closer 
to data thanks to a shift to lower px- For example for the V2 of pions gets closer to a 
scaling with mx — mo as found experimentally [70| , which can be seen as an indirect effect 
of energy conservation. More recently energy conservation in the coalescence process has 
been investigated by Ravagli and Rapp 69J using a formulation of coalescence from the 
Boltzmann collision integral. They find that if energy conservation is implemented then the 
V2 scaling as a function of px moves toward a scaling versus mx — mo . Another important 
development is the inclusion of a mass distribution for the quarks as an effective way for 
including in-medium interaction i6j, this is another way to allow for both momentum and 
energy conservation and leads to a fairly good agreement with the data, for more details see 
the contribution of T. Biro to these Proceedings. 

An issue related also to energy conservation is the entropy conservation. Coalescence reduces 
the number of particles by about a factor two (if resonance production is not taken into 
account), this rises suspects of a strong decrease of entropy. Of course entropy is not the 
number of particle, but depends also on the degeneracy of the two phases and on the mass 
of the particles. For example converting one gluon into one pion with the same mass leads 
to a huge entropy violation. In other words a large chemical potential is enforced by number 
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^ The V2{pt) used here is only schematic to show the effect of the violation that for a more realistic 
V2{pt) is however reduced to about a 10% for Ap — 0.24 GeV. 
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conservation causing an entropy decrease unless at the same time there is a large volume 
expansion. 

In Fig.[5]it is shown the evolution of dS/dy (for \y\ < 0.5 as a function of time if one assumes 
that coalescence takes place during a mixed phase of 4 fm (and 8 fm upper line), during which 
the volume expands with P = 0.5 and quarks are gradually converted into hadrons which are 
allowed to decay according to their width. For the full calculation (expansion plus decays) a 
decrease of about 15 % is found, but in our coalescence model there is also a decrease in the 
energy by about 15% {S = {E + lnZ)/T, where Z is the partition function). Therefore if 
coalescence is implemented taking care of the energy conservation probably entropy can also 
be conserved. For the multiplicity of various hadrons the calculation is based on the results 
from our coalescence model (but only u,d matter is considered). The important role of the 
decays and the moderate role of the expansion are shown by the different cases considered, 
see Fig. [5l We also see that assuming a longer mixed phase (i/ = 8 fm), and therefore 
a larger volume at the end entropy can be conserved even in such a simple approach. In 
addition one should also take into account the interaction among quarks; it has been shown 
using the IQCD equation of state for an isentropically expanding fireball that the evolution 
of the effective number of particles reduces significantly around the crossover temperature 
|71j . This of course helps to solve the entropy problem inherent to quark coalescence, as 
pointed out also by Nonaka et al. [51] . 

4. Space-momentum correlation and V2 scaling - The constituent quark number scaling of 
the elliptic flow was derived by Voloshin and Molnar assuming that the x-space can be 
decoupled from p-space and integrated out. This means that the scaling has been explicitly 
demonstrated only if the coalescence probability is omogeneous in space and if space and 
momentum are not correlated. A detailed discussion on effects coming from space and 
momentum correlation as a source of scaling violation can be found in Ref . j46| together 
with classes of distribution that lead to an approximate scaling. Another study on the effect 
of phase space distribution can be found in Ref. [66j . 

Results presented here with our coalescence model do not assume the coalescence scaling 
and show that if radial flow space-momentum correlation (but still in a uniform fireball) is 
included the scaling is still recovered [S]. An important tool of investigation in this context 
is supplied by parton cascade studies that calculate the time evolution of the phase space. 
A first investigation by Molnar finds that the scaling between baryons and mesons still 
persists even if a strong violation is found respect to the V2 at quark level [S3]. However it 
is not clear the dependence on the freeze-out criteria, on the wave function width and on 
the interplay with indipendent fragmentation mechanism as sources of scaling breaking, see 
also D. Molnar contribution to these Proceedings. 

Finally it is clear that a better understanding of HBT measurements can supply fundamental 
information on this issue that is the potential main source of scaling violation and therefore 
should be investigated in deeper detail. At LHC the different dynamical evolution may lead 
to an even stronger scaling violation. 

5. Jet-like correlation - The last but not the least issue is related to the correlation among 
hadrons in the hadronization process. At variance with the other issues this is the one 
more relavant at intermediate pr- At RHIC it has been possible to measure the correlation 
between hadrons at momentum p^*^ >4-6 GeV and the associated particles with momenta 
Pt < Px^^ as a function of the azimuthal angle respect to the trigger particle |75|76] . Such 
a kind of measurement has shown that hadrons at intermediate pr come with associated 
particles at lower momentum in the opposite direction. In the seminal papers |8|9j we have 
taken into account the possibility of coalescence between thermal partons and minijets and 
this was a primitive way to propose the idea that at intermediate pr particle from coalescence 
could follow the di-jet correlation. A formal framework for the effect of quark correlations 
has been studied by R. Fries et al. [57] and R. Hwa and C.B. Yang (62j. They have shown 
that coalescence goes along with particle angular correlation, even if more stringent tests 
are currently under investigation. Anyway the main problem remains the study of the origin 
of correlation at parton level which asks for a transport approach, a first explorative study 
has been performed by Molnar [64] . 
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In the future it is desireable to have a deeper theoretical investigation of coalescence, in fact 
it is easy to predict that such a kind of processes will be even more dominant at the Heavy-Ion 
experiments at the Large-Hadron-Collider (LHC) at CERN, see also P.Levai contribution to 
these Proceedings. In fact hard processes will be largely dominant together with a larger jet 
quenching and a thermal source with larger radial flow which should lead to an overwhelming 
production of hadrons from coalescence in a wider momentum range |79l80j . Therefore the 
features of more sophisticated models will have wide chances to undergo many experimental 
tests at both the LHC and RHIC programs. 



6 Phase Space Coalescence for heavy quarks 

The coalescence model has been also applied to study open and hidden charm and bottom 
production. For the evaluation of charmonia yield the idea of coalescence j58j or recombination 
at phase boundary 59] has been investigated by several authors as a competing effect to char- 
monia suppression. We have suggested instead that the dominance of coalescence mechanism in 
a wide range of momenta can be applied also for the study of the D and B mesons |10|17|89|90j 
that are the principal tool for investigating heavy quark interaction and thermalization in the 
QGP. In a first paper we have shown that such an approach allows to estimate both spectra 
and elliptic flow of the D mesons and that this offers the possibility to relate J/tf' and D meson 
spectra if both come mainly from coalescence [10] . 

Surprisingly, data from the Relativistic Heavy-Ion Collider (RHIC) for single electrons (e^) 
associated with semileptonic B and D decays in semi-central Au-Au collisions exhibited a V2 of 
up to 10% [82|85j . indicating substantial collective behavior of charm (c) quarks consistent with 
the assumption of a c-quark V2 similar to the one for light quarks, apart from a px shift due to 
radial- flow effects [TU] . In addition, the nuclear suppression factor was found to be comparable 
to the pion one, 0.3 |83|84j . 

On the other hand for heavy quarks (HQ) the predicted amount of energy-loss has been 
shown to be insufficient (at variance with the light quark case) to account for the observed non- 
photonic single electrons nuclear suppression (small Raa) and for a strong degree of collectivity 
(large 1)2) [5S^. Therefore for HQ the challenge is mainly the understanding of the quark in- 
medium interaction, even if the acquired knoweledge on hadronization mechanism from light 
quarks plays a significant role. 

Based on lattice QCD (IQCD) which suggests a resonance structure in the meson-correlation 
function at moderate temperatures |77| , an effective model for heavy- light quark scattering via 
D and B resonances has been suggested. In Ref.[16] it was found that resonances cause a 
reduction of about a factor of three in the thermalization time respect to pQCD estimates. 
To relate such a microscopic description to the single elctron data, we employ a Fokker-Planck 
approach for c and b quarks in the QGP based on elastic scattering with light quarks via D- and 
_B-meson resonances with a width _r=400-750 MeV (supplemented by perturbative interactions 
in color singlet channel) [16]. Heavy-quark (HQ) kinetics in the QGP is then treated as a 
relativistic Langevin process |17| : 

dfg _ d{pfQ) d^fQ 

where /q is the HQ phase space distribution, 7 and Dp are the corresponding drag and (momen- 
tum) diffusion constants which determine the approach to equilibrium and satisfy the Einstein 
relation, T ~ Dp/^Mq. The medium is modeled by a spatially homogeneous elliptic thermal 
fireball which expands isentropically. Finally the hadronization is treated by a coalescence (see 
Eq.([T])) plus fragmentation approach with the distribution of HQ that undergoes a fragmen- 
tation process evaluated as /cApt) * [1 - Pc,b->{D,A,),(B^At){PT)], where Pc,b->{D,A,)XB,At) is 
the probability for a HQ to coalesce according to Eq.jT]). Results for Au-|-Au at 200 AGeV 
from the Langevin simulation including hadronization by coalescence-f fragmentation (left) and 
fragmentation only (right) are shown in Fig|6] together with experimental data [83186] . It is 
clear that elastic scattering in a pQCD scheme is insufficient to account for the small Raa and 
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Fig. 6. Nuclear suppression factor for single-electrons, including coalescence and fragmentation at 
hadronization (left panel) and only with fragmentation (right panel), see text. 



large V2 (see FigjT]) independently on the hadronization scheme applied. The red band shows 
the full calculation with c, b quarks that scatter in the presence of hadronic-like resonances with 
a width 7^ ~ 0.4 — 0.75 GeV. We notice that contamination of single-electrons from B decay 
is significant already at pt ~ 2 GeV (corresponding to a cross-point between c and b spectra 
around 4—5 GeV in agreement with FNNLO calculation [57]). Therefore it is necessary to in- 
clude the B mesons (despite the inherent uncertainties in the b/c ratio) to draw any conclusion 
on the interaction processes behind the experimental results. 
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Fig. 7. Elliptic flow for single-electrons, including coalescence and fragmentation at hadronization (left 
panel) and only with fragmentation (right panel), see text. 



Comparing the band in the left and right panel of Figs. [6] and [71 we can appreciate the effect 
of coalescence in the increasing of Raa at ~ 1 — 4 GeV and the simultaneous increase of the 
elliptic flow V2, see FiglTl (9110] . Therefore coalescence mechanism reverts the usual correlations 
between Raa and V2 , and allows for a reasoble agreement with the experimental data of both 
Raa and V2. 

Very recently within a Brueckner many-body scheme in medium T-matrices for heavy-quarks 
scattering off light quarks |89J have been evaluated starting from IQCD potential. Therefore 
the existence and width of D and i?-like resonances are no longer assumed like in the effective 
lagrangian but extracted from IQCD. It is found even a better agreement with the data [55] . 
Even if inherent uncertainties in the extrapolation of the potential have to be evaluated in the 
next future, this approach constitutes a promising tool to connect the observables to the infor- 
mation from IQCD. Furthermore we note that if resonant scattering with increasing strength at 
decreasing temperature is the dominat interaction channel this would lead to a natural merg- 
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ing into a quark-coalescence. Therefore heavy quarks may aUow a coherent description between 
in-medium interaction that drives the thermahzation process and the subsequent hadronization. 

We also mention that LHC will play an important role to validate such a picture. In Ref. 
[50] predictions for LHC have presented. They are quantitatively rather similar to our RHIC 
results |17j , due to a combination of harder initial HQ-pT spectra and a decrease in interaction 
strength in the early phases where non-perturbative resonance scattering is inoperative. There- 
fore we conclude that if at RHIC the dominant contribution to HQ interactions are hadron-like 
resonances, at LHC we should observe an Raa and V2 pattern similar to RHIC. 

Finally we point out that in the next future we can look not only at the quarkonia yields, but 
scrutinize their pT which allows to check the self-consistency between D and B mesons and J/'F 
and T quarkonia that in the QGP should be related by the unique underlying HQ distribution. 
Therefore phase space coalescence will provide a much deeper insight into the long-standing 
issue of J/^ suppression and regeneration. 



7 Summary and Conclusions 

The first stage of RHIC program has shown clear signs of modification of the hadronization 
mechanism respect to pp collisions in the light quark sector. There are several evidences that 
hadronization proceeds through coalescence of massive quarks close in phase space. The baryon- 
to-meson ratio enhancement observed at intermediate pt and the scaling of the elliptic flow 
with the number of valence quarks, are robust effects against improvements of the naive coa- 
lescence picture: inclusion of resonances, finite width of the wave function, effective quark mass 
distribution, gluons in the higher Fock staes of the wave function, energy conservation ... All 
these studies that improve the original coalescence picture from a theoretical point of view 
lead generally to an even better agreeement with the data. The main open question remains 
the development of a dynamical coalescence model that can shed light on the issue of space- 
momentum correlation and jet-like correlations that could be significant sources of deviation 
from naive coalescence. 

The knoweledge of hadronization for light quarks seems to play a role also in the new chal- 
lenges posed by heavy-quark probes. Here the main issue is the dominant interaction mechanism 
and its relation to the microscopic structure of the QGP. The role of coalescence is mainly a 
modification of the correlation between Raa and V2. Finally, coalescence in the heavy quark 
sector can provide an inherent consistency between the in-medium interaction of HQ and the 
subsequent hadronization considering that a pole in the HQ propagator above Tc can be viewed 
as a precursor of coalescence. 
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